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Abstract 



We present a next-to-leading order QCD calculation of the cross section and 
longitudinal spin asymmetry in single-inclusive charged-lepton production, pp — )■ 
£^X, at RHIC, where the lepton is produced in the decay of an electroweak 
gauge boson. Our calculation is presented in terms of a multi-purpose Monte- 
Carlo integration program that may be readily used to include experimental spin 
asymmetry data in a global analysis of helicity parton densities. We perform a toy 
global analysis, studying the impact of anticipated RHIC data on our knowledge 
about the polarized anti-quark distributions. 



1 Introduction 



Despite much progress over the past three decades, many open questions concerning the hehcity 
structure of the nucleon remain. For example, we so far have only a rather sketchy picture 
of the individual polarizations of the light quarks and anti-quarks, Am/m, Am/m, Arf/rf, Aj/d, 
where the helicity parton distributions are as usual denoted by Ag, Ag, and their spin-averaged 
counterparts by g, q. Since nucleons have up and down quarks as their valence quarks, the 
light-quark and anti-quark polarizations are of much interest in QCD and play key roles in 
many models of nucleon structure and, more generally, for our fundamental understanding of 
the nucleon pQ. 

While lepton scattering has provided fairly precise and solid information on the "total" up 
and down distributions. Am + Am, Ad + Arf, through inclusive measurements, information on 
the individual Am, Am, Ad, Ad is much more sparse and still afflicted by large uncertainties |2l 
|3]. The tool exploited here so far has been semi-inclusive deep- inelastic scattering (SIDIS), 
in which one detects a specific hadron in the final state and uses the fact that the fiavor 
content of that hadron will typically be correlated with the fiavor of the quark or anti-quark 
hit by the virtual photon in the basic deep-inelastic reaction. Measurements of SIDIS spin 
asymmetries have vastly improved in recent years. By now, quite precise data sets are available 
for various differently produced hadrons [IHS]. On the other hand, extraction of polarized 
parton distributions from SIDIS relies on the applicability of a leading-twist factorized QCD 
description of the reaction, allowing among other things the use of fragmentation functions [9] 
for the produced hadron determined from other processes. For the kinematics accessible in 
SIDIS so far one may worry if, for example, sub-leading twist effects can really be ignored in 
the theoretical analysis of the data. This leads to an uncertainty that is presently hard to 
quantify. It does have to be said, however, that probably for the light quarks and anti-quarks, 
which are primarily determined from pion SIDIS, this is less of a concern. Also, thanks to the 
recent COMPASS measurements the kinematic reach of SIDIS data has become quite large now, 
extending into a regime where higher-twist effects should be less relevant. In any case, as with 
any measurements of nucleon structure, it is of great value to have a completely independent 
probe that does not involve any hadronic fragmentation and is characterized by momentum 
scales so large that perturbative calculations are expected to be reliable and sub-leading twist 
effects unimportant. 

It has long been recognized that boson production at the Relativistic Heavy Ion Collider 
RHIC may provide unique and clean access to the individual helicity polarizations of quarks 
and anti-quarks in the colliding protons Thanks to maximal violation of parity in the 

elementary Wqq' vertex, W bosons couple to left-handed quarks and right-handed anti-quarks 
and hence offer direct probes of their respective helicity distributions in the nucleon. Since 
spin asymmetries obtained from a single longitudinally polarized proton beam colliding with 
an unpolarized beam are parity-violating for sufficiently inclusive processes, they have become 
the prime observables in the physics program with W bosons at RHIC [T014T2] : 
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Here the a~^^ etc. denote cross sections for scattering of protons with definite hehcities as 
indicated by the superscripts. As one can see, the hehcities of the second proton are summed 
over, leading to the single-spin process pp — )■ W^X. The basic idea behind measurements of 
the helicity distributions at RHIC is then as follows: production of W~ , for example, selects 
primarily a u anti-quark from one proton in conjunction with a (i-quark from the other. Thus, 
for the simple lowest-order (LO) parton-model process du — ?• W~ the single-spin asymmetry 
becomes 

where Ac?, Am denote the usual helicity distributions, probed here at a scale of the order of the 
W mass Mw, and where 

42 = ^e*- (3) 



with the rapdity yw of the W bosons and the hadronic center-of-mass energy V S. It follows 
that at large y, where x'J ~ 1 and <^ 1, the asymmetry will be dominated by the valence 
distribution probed at x^, and give direct access to —Ad{xi,M^)/d{xi,M^). Likewise, for 
large negative y, is given by Au{x^)/u{x^). The situation for follows analogously. 

In practice, the above reasoning needs to be augmented in various ways. Foremost, there 
is an experimental issue: the detectors at RHIC are not hermetic, which means that missing- 
momentum techniques for the charged-lepton plus neutrino (iue) final states cannot be used to 
detect the W and reconstruct its momentum. Instead, the strategy adopted by the RHIC ex- 
periments is to detect the charged decay lepton and determine its transverse momentum p^i and 
rapidity rji. The relevant process therefore becomes the single- inclusive reaction pp — )■ i^X, 
similar in spirit to the processes pp — )■ tcX, pp — )■ jetX [T51 - IT7] used at RHIC to determine 
gluon polarization in the nucleon. The ensuing expression for the single-spin asymmetry for 
pp £^X differs from that in Eq. (|2]) even at lowest order, since the lepton transverse momen- 
tum and rapidity do not completely determine the momentum fractions of the initial partons, so 
that an integration over momentum fractions appears in the expression. Nevertheless, studies 
have shown [T2Hll] that despite this fact there should still be excellent sensitivity to the distri- 
butions Am, Am, Ad, Ad for appropriately chosen lepton kinematics. Very recently, the RHIC 
collaborations have presented the first preliminary data on the cross section and single-spin 
asymmetry for boson production at RHIC [T8|[T9]. 

There are also theoretical issues that modify the simple picture given by Eq. ([2]), regardless 
of whether one uses W or lepton kinematics. Cabibbo-suppressed contributions, which involve 
the polarized and unpolarized strange quark distributions, and also contributions by Z bosons, 
are relatively straightforward to take into account. A more important issue is the higher-order 
QCD corrections to the leading-order (LO) process qq' — )• W^. At next-to- leading order (NLO), 
for example, one has the partonic reactions qq' — )■ VV^g and qg — )■ W^q'. Despite the fact that 
the W mass sets a rather large scale so that the strong coupling constant as{M]y) is not large, 
the corrections can be significant and certainly need to be known for a reliable theoretical 
extraction of spin-dependent parton distributions from data. Consequently, a lot of theoretical 
work has gone into the calculation of higher-order QCD corrections to the spin asymmetries in 
weak-boson production at RHIC. Early work in this area |20] - [26] focused on the case where the 
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W boson is detected directly, which is kinematically simpler and allows one to obtain analytical 
results for the NLO corrections. More recently, also an all-order resummation of terms in the 
partonic cross section that are logarithmically enhanced near partonic threshold was presented 
for this case [27]. While, as we discussed above, the direct detection of the W kinematics is not 
possible at RHIC, the relative size of the NLO corrections is expected to be rather independent 
of whether one takes into account the W — >■ iu^ decay or not, since this decay does not involve 
any strong interactions and all QCD corrections occur only in the initial partonic state. 

There have also been extensive studies of higher-order QCD corrections for the experimen- 
tally more relevant case of single- inclusive lepton production, pp — )■ i^X. In Refs. [13 1 [28] 
the program RHICBOS was introduced. RHICBOS is a Monte-Carlo integration program for 
lepton distributions, specifically adapted to the polarized pp collisions at RHIC. It puts par- 
ticular emphasis on the effects of soft-gluon emission and their impact on the region when the 
produced intermediate vector boson has small transverse momentum, qt- In the lowest-order 
diagrams qq' ^ W ^ iue, one has gr = 0. Gluon radiation generates a recoil transverse mo- 
mentum. When qT tends to zero, large logarithmic corrections develop in the spectrum of the 
W^s. These can be resummed to all orders in perturbation theory, following the Collins-Soper- 
Sterman (CSS) formalism |22j, which is done in RHICBOS at next-to- leading logarithmic level. 
RHICBOS is widely used for phenomenological studies related to the W program at RHIC (see, 
for example, Ref. [12]). 

Despite this earlier work, we will present in the present paper a new NLO calculation of the 
cross sections and spin asymmetries for pp — )■ i^X at RHIC. There are several reasons why this 
is in our view a necessary addition. First of all, it is of value to have an independent calculation 
of the relevant observables at RHIC. Second, the RHIC data for the spin asymmetry in W 
production will ultimately need to be included in a global NLO analysis of parton distributions 
that includes all available information from lepton scattering and pp collisions at RHIC. Only 
then can the best possible information on the Aq and Aq be extracted. Inclusion of pp scattering 
data in a global analysis is a relatively complex task [2j since the computation of the parton 
subprocess cross sections is typically numerically quite involved. Recent papers [21[3| used a 
technique based on Mellin moments [30] to achieve the first global analysis of polarized lepton 
scattering data and data for pp — )■ ttX and pp — )■ jetX from RHIC. Our calculation presented 
in this paper is set up in such a way that inclusion of RHIC data for pp — )■ i^X will be 
straightforward. This is an advantage over RHICBOS, which requires prior computations of 
certain "grid" files for a given set of parton distribution functions, and hence is to our knowledge 
not readily suited for use in a global analysis code. Using our new code, we will present in this 
paper a first "toy" study of a global analysis that includes projected or estimated data for W 
observables at RHIC. 

Finally, we also have a more theory-related reason for performing a new NLO computation 
of the W observables at RHIC. As described above, RHICBOS includes the all-order resum- 
mation of large logarithmic corrections arising at small W transverse momentum q^- As is 
well-known [31j, these logarithms are very relevant if one is interested, for example, in the 
low-transverse momentum distribution of the W boson itself. However, for the single-inclusive 
lepton cross section, the situation is somewhat different, g^-- resummation is really only useful 
when the observable is directly sensitive to (small) q^- For pp — )■ i^X this is the case when the 
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measured lepton transverse momentum, is in the vicinity of Mvy/2. This may be understood 
as follows: for the LO reaction qq' ^ W ^ iui, the lepton transverse momenta are basically 
limited to pj; < Mwf^, except for effects related to the finite decay width of the W. This 
means that lepton transverse momenta > M\y/2 primarily arise from higher-order gluon 
emission. Just above the lepton transverse momentum spectrum is then dominated by 

the same logarithms that are present in the W transverse momentum distribution, which are 
resummed by the CSS formalism. This is the motivation behind the resummation implemented 
in RHICBOS. In practice, however, the RHIC experiments sample over a fairly broad range 
of ■ For the theoretical calculation this implies integration of the observables over pj; over 
this range. The broader this range, the less dominant are the soft-gluon effects addressed by 
resummation, and the less useful is resummation. This becomes particularly evident for the 
rapidity-distribution of the lepton, integrated over p^i , which is the most relevant observable in 
VT-physics at RHIC and, in the spirit of Eq. ([2]), the best tool to separate the various polarized 
parton distributions. To state our point more succinctly: for practical purposes at RHIC, the 
region j9j; ^ M\y/2 is only a relatively small part of the sampled kinematics, so that the qx 
logarithms are not expected to dominate the observables. Their resummation is, then, not re- 
ally appropriate and does not necessarily lead to an improvement of the theoretical calculation, 
since at the level of pr, -integrated observables there will be other higher-order effects that are 
of the same size as those provided by the terms logarithmic in qx [331435] . In any case, from 
the point of view of extracting polarized parton distributions, it is advisable in our view to 
select observables at RHIC that are insensitive to the complications associated with multiple 
soft-gluon emission. It therefore seems preferable to us to use a plain NLO calculation for 
studies of VT-production at RHIC, which we aim to do in this paper. 

The remainder of this paper is organized as follows: in the next section we very briefly 
discuss our NLO calculation, which is overall quite standard. The main part of the paper is 
then phenomenological and presented in the following two sections. Apart from providing NLO 
predictions for RHIC in Sec. [31 we also present in Sec. [Ha "proof-of-principle" study of a global 
analysis involving W spin asymmetries. We finally conclude in Sec. [5] 



2 Next-to-leading Order Calculation 

In order to evaluate the NLO QCD corrections to the process we rely on the version of the 
subtraction method introduced and extensively discussed in Refs. [MIISZ], and later extended 
to the polarized case in Ref. [16]. We refer the reader to those references for the details. 
Figure [1] shows some of the Feynman diagrams contributing at LO and NLO in the case of W 
exchange. The calculation is implemented in the Monte-Carlo like code 'CHE' (standing for 
'Collisions at High Energies') which provides access to the full kinematics of the final-state 
particles, allowing the computation of any infrared-safe observable in hadronic collisions and 
the implementation of realistic experimental cuts. It is worth noticing that the same code can 
compute the unpolarized, the single polarized and the double polarized cross sections. 

^The code is available upon request from deflo@df.uba.ar 
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Figure 1: Feynman diagrams for W production with leptonic decay: (a) leading- order, (b) NLO 
virtual correction, (c) NLO real emission, (d) NLO quark-gluon scattering. Crossed diagrams 
are not shown. 

Besides the contribution driven by W exchange, the code also allows the computation of the 
background arising from Z-boson and/or photon exchange at the same accuracy in perturbative 
QCD. We point out that at NLO the contribution from photon exchange, qq — )■ '~f*g followed 
by 7* — )■ , may generate large contributions when the high-transverse momentum photon 
splits almost collinearly into the lepton pair, producing high-pT; leptons with a very low invariant 
mass. A proper treatment of this configuration would require the addition of a fragmentation 
contribution based on parton-to-dilepton fragmentation functions |38] . On the other hand, it is 
likely that configurations with two nearly coUinear leptons can be distinguished experimentally 
from true single high-pj; leptons. In our calculation we avoid such configurations by requiring 
the lepton pair to have an invariant mass Mi^^i^ > 10 GeV. 

We note that we have checked the results for the spin-averaged cross section in our code 
against the MCFM [39] and DYNNLO [lO] codes. We have also computed the fully inclusive 
spin-averaged and polarized W cross sections, integrated over all lepton angles. For these cross 
section, analytical results are available [231425] . with which we agree. 

A virtue of our code is that it lends itself to inclusion in a global analysis of polarized 
parton distributions, along the lines presented in Refs. [2l|3]. In these papers, a method based 
on Mellin moments was used [30], for which the theoretical expression for any cross section 
is split up into parts that are independent of the parton distributions, coupled to the Mellin 
moments of the parton distributions. This procedure was shown to tremendously speed up 
the NLO fit, since the pieces that do not depend on the parton distributions, which contain 
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the most time-consuming computations, can be calculated "once and for all" prior to the fit 
and stored as large arrays. In the actual fit one then only needs to perform numerical inverse 
Mellin transforms, which is straightforward. As was shown in Ref. [2], the computation of 
the pre-calculated factors is possible in a timely manner even for a code based on Monte- 
Carlo integration, provided a proper importance-sampling is used. We have implemented the 
corresponding strategies described in [2] in our code. 



3 Phenomenological Results for RHIC 

We now use our NLO code to present some numerical results for polarized pp collisions at 
RHIC at center-of-mass energy = 500 GeV. Our default choice for the spin- dependent 
parton distribution functions is the DSSV set [2l|3]. Since we want to study the sensitivity of 
different observables to the polarized parton distributions, we will also consider a few other 
(and less recent) sets of polarized densities that primarily differ in the anti-quark polarizations: 
the "standard" and "valence" sets from GRSV [41J, which have SU(2) symmetric and broken 
sea distributions, respectively, and the " DNS-Kretzer" and "DNS-KKP" sets |12]. The last two 
sets correspond to fits to the same data for inclusive and semi-inclusive lepton scattering, but 
obtained using different sets of fragmentation functions [ISJIH] to analyze the semi- inclusive 
asymmetries. We note that not all of these additional sets of polarized parton distributions 
are completely compatible with all information now available from SIDIS. However, given the 
potential uncertainties in SIDIS mentioned in the Introduction, they are useful in order to 
gauge the sensitivity of future RHIC measurements. Figure |2] shows a comparison of the light 
anti-quark distributions xAu{x, Q) and xAd{x, Q) for the various sets, evaluated at the scale 
Q = 80 GeV relevant for W production. As one can see, there are large differences among 
the distributions, both in qualitative behavior regarding breaking of flavor-SU(2) symmetry, 
and in magnitude. Since all sets provide very similar results for the sum of quark and anti- 
quark polarized distributions, Agj + Ag^, and since the observables are rather insensitive to the 
polarized gluon density, differences in the W spin asymmetries computed with the different sets 
can be mostly attributed to the differences in the sea distributions. 

In order to compute the unpolarized cross section in the denominator of the spin asym- 
metries, we use the MRST2002 NLO set. This choice is motivated by the fact that 
the MRST2002 set was also used as the "baseline" unpolarized set for the DSSV spin- 
dependent parton distributions. We have verified that the use of more recent sets of distribution 
functions, like CTEQ6 |16], results in very similar cross sections. 

We set the masses of the vector bosons to Mz = 91.1876 GeV and Mw = 80.398 GeV, 
and the corresponding decay widths to Tz = 2.4952 GeV and Tw = 2.141 GeV [S]. We 
neglect contributions from b and t quark initial states to W production and, consistent with 
that, use the following values for the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements: 
\Vud\ = \Vcs\ = 0.975 and \Vus\ = \Vcd\ = 0.222. We do not include any QED or electroweak 
corrections, but choose the coupling constants a and sin^ 9w in the spirit of the "improved Born 
approximation" |l8lll9], in order to effectively take into account the electroweak corrections. 
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Figure 2: Left: Next-to-leading order xAu{x,Q) evaluated at the scale Q = 80 GeV for the 
DSSV IMW (solid), GRSV gT]/ (dashed), DNS-Kretzer (dotted) and DNS-KKP (short- 
dashed) sets of polarized pdfs. Right: Same for xAd{x,Q) (right-hand side) 



This approach results in sin'^ 9w = 0.23119 and a = a{Mz) = 1/128. We also require the 
lepton pair to have an invariant mass Mi-^i,^ > 10 GeV, in order to avoid potentially large NLO 
contributions from production of a high-p^ nearly real photon that subsequently decays into a 
pair of almost collinear leptons, as discussed in Sec. O 

We will study two different observables for lepton production in pp — )■ i^X: the transverse 
momentum (pj,) distribution of the charged lepton with a rapidity cut of \rii\ < 1, and the 
rapidity distribution with pxi > 20 GeV. We count rapidity as positive in the forward direction 
of the polarized proton. There are two hard scales in the process, which are of the same order: 
the mass of the gauge boson and the transverse momentum of the observed lepton. We choose 

= = (M^ +p|.J/4 as the default factorization and renormalization scales for both the 
W and the Z/^ contributions. We note that the scale dependence of the cross sections and, in 
particular, the spin asymmetries is extremely mild in case of vector boson production, so that 
other choices like fip = fJ'R = Mw give very similar results. 

We start by investigating the contribution by Z and 7 exchange to the cross-sections. 
Even though the qqZ coupling is also parity-violating and hence may contribute to the single- 
longitudinal spin asymmetry, the Z and 7 contributions are to be regarded in a sense as "dilu- 
tions" of the W signal. Being almost symmetric in and i~ , they will somewhat decrease the 
clear-cut sensitivity of Al to the polarized sea-quark distributions, and they also contribute to 
the spin-averaged cross section. Figure [3] compares the Z/ 7 contributions (dashes) to the ones 
arising from W (solid), for the spin-averaged cross section as functions of rapidity (left-hand 
side) and transverse momentum (right-hand side). As can be observed, for positrons (upper 
row) the contribution from Z turns out to be rather small in the central rapidity range (about 
~ 7%) while it does significantly add for electrons (lower row), reaching more than 40% of the 
W~ contribution. The effect turns out to be more noticeable at larger rapidities. We note that 
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Figure 3: Contribution from W (solid lines) and \Z\^ + 7Z + I7P (dashed lines) production 
to the rapidity (left-hand side) and transverse momentum (right-hand side) distributions of the 
leptons in spin- averaged collisions at RHIC. The upper plot corresponds to positron and the 
lower plot to electron production in unpolarized collisions. 



the |Zp contribution dominates strongly over the 7Z interference and I7P ones. In the case 
of the transverse momentum distribution, the relative contribution by Z's strongly depends on 
PTi- Since the peak of the distribution occurs around pti ~ Mv/2, with My the vector boson 
mass, the difference between the W and Z masses induces dominance of the Z contribution for 
PT, > 45 GeV. 

The same comparison is shown in Fig. IHfor the single-polarized case, where we rely on the 
DSSV set of parton distributions for the polarized beam. While Z exchange can produce a 
single-spin cross-section, its parity violation component is rather small and results in a contri- 
bution that does not exceed a few percent of the dominant W one. Of course, this does not 
mean that the Z contribution can be neglected in the analysis to be carried out to extract the 
polarized parton distributions at RHIC. The main observable is the spin asymmetry which, at 
least in the case of electron production, can be considerably reduced by the Z/j contribution in 
the unpolarized cross section. Vetoing events with a lepton pair could be valuable in this case 
in order to increase the sensitivity to the polarized parton densities. Conversely, Z bosons by 
themselves may offer an interesting advantage if both charged leptons from the decay Z — > 
can be detected, because in that case one can in principle reconstruct the kinematics of the Z 
boson, which is not possible for the Ws because of the neutrino in its decay. One would then 
be able to directly access the momentum fractions of the parton distributions, in the spirit of 
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Figure 4: Same as Figl^ but for single-polarized collisions using the DSSV parton distributions. 



Eq. ([2]). Unfortunately, statistics for reconstructed Z decays with both decay leptons will likely 
remain rather low at RHIC. 

We next investigate the kinematics of W production at RHIC. In our view, it is preferable 
to consider distributions in lepton rapidity, rather than transverse momentum, since there is a 
particularly strong and direct correlation between lepton rapidity and the partonic momentum 
fractions. This correlation was already evident in the LO asymmetry as a function of the W^s 
rapidity discussed in Eq. (|2]), for which we had Xi_2 = Mm^^vw _ Qne can expect that, at least 
to some extent, this relation between momentum fractions and rapidity at the gauge boson 
level will be inherited by the lepton. Figure E] shows the correlation between the averages of the 
momentum fractions, {xi^2), and the rapidity of the charged lepton computed at NLO accuracy 
for W~ (left-hand side) and production (right-hand side) in spin-averaged collision^. A 
remarkably strong correlation is found between (xi,2) and rji in both cases. Large negative 
lepton rapidity corresponds to small (large) momentum fractions Xi (2:2). The opposite occurs 
for large positive rapidities. Actually, as a rough approximation one can parameterize these 
correlations by the simple "empirical" formulas 



(4) 



Considering that RHIC experiments will allow one to reach rapidities of the order of |?7i| ~ 2, 



^The correlation remains the same when the Z/7 contribution is included. 
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one can expect sensitivity to the polarized quark and anti-quark distributions in the region 
0.05 < a; < 0.4. We note that similar results as in Fig. [5] were also found in Ref. [H]. 




Figure 5: Averages of the momentum fractions xi^2 as functions of the charged lepton's rapidity 
Tji for W~ (left) and production (right) at RHIC. 

Because of the correlation shown in Fig. [5l the combinations of parton distributions pre- 
dominantly probed will vary with rji. However, here also the underlying structure of the weak 
interactions enters. For W~ production, neglecting all partonic processes but the dominant 
ud — )■ W~ — )• e~z/e one, the asymmetry is found to be given by 



^e- _ 4(xi,X2) [^u{Xi)d{x2){l - COsOf - M{xi)u{x2){l + COS (, , , 

^ ^ [^(^i)^(^2)(l - cos^)2 + d(a;i)M(x2)(l + COS^)2] 



where /^(^^ denotes an appropriate convolution over momentum fractions, and where 6 is 
the polar angle of the electron in the partonic c.m.s., with ^ > in the forward direction 
of the polarized parton. Note that 6 itself depends on the momentum fractions and on the 
lepton's rapidity. At large negative rji, one has X2 ^ xi and 6* ~ tt. In this case, the first 
terms in the numerator and denominator of Eq. ([5]) strongly dominate, since the combination 
of parton distributions, Au{xi)d{x2), and the angular factor, (1 — cos 6*)^, each dominate over 
their counterpart in the second term. Therefore, the asymmetry provides a clean probe of 
Au{xi)/u{xi) at medium values of Xi. By similar reasoning, at forward rapidity rji ^ 
the second terms in the numerator and denominator of Eq. ([5]) dominate, giving access to 
—Ad{xi)/d{xi) at relatively high xi. For the production channel one has instead of ([5]) 

_ [Ad{x,)u{x2){l + cosOy - An(xi)rf>2)(l - cosOy] 

^ /®(xi,x2) [d{xi)u{x2){l + cosOy + u{xi)d{x2){l - cosOy] 

Here the distinction of the two contributions by considering large negative or positive lepton 
rapidities is less clear-cut than in the case of W~ . For example, at negative r]i the partonic 
combination d{xi)u{x2) will dominate, but at the same time 6' ~ vr so that the angular factor 
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(1 + COS 6')^ given by the basic electroweak interaction is small. Likewise, at positive t]i the 
dominant partonic combination Au{xi)d{x2) is suppressed by the angular factor because 6* ~ 0. 
So both terms in Eq. (Q will compete essentially for all rji of interest. This is also the reason 
why the W~ cross section can become larger than the one at high rapidities (see Fig. E]). 
As was discussed in Refs. [SUIEI], the study of hadronic W decays, in particular of charmed 
final states [20], might help in improving this situation. 



The features displayed by Eqs. (IS]),© are fully reflected in the behavior of the calculated 
NLO spin asymmetries. Figure |6] shows Aj^ for electrons and A'l for positrons at RHIC, 
as functions of the charged lepton's rapidity, considering only leptons arising from boson 
exchange. We are now using all of the sets of polarized parton distributions that we introduced 
earlier. The spread in the predictions for the asymmetry at r/; < directly reflects the 
dispersion in both the absolute magnitude and sign of the different Au{x) distributions shown 
in Fig. [2] in the range 0.05 ^ a; < 0.2. On the other hand, the asymmetry becomes large and 
positive at high rji, which reflects the fact that Ad{x) remains negative at high x for all sets 
of polarized parton distributions considered here. A'l does not show as clear features, for the 
reasons we just discussed. Nevertheless, at r/; > one can observe again that the spread of 
the predictions for the asymmetry is quite strongly correlated to the one found for the Ad 
distributions at 0.15 ^ x < 0.3 in Fig. [2l Overall, the asymmetry is negative because of the 
contribution from Am in Eq. (Q, which is known to be positive from lepton scattering. 
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Figure 6: Rapidity dependence of the NLO single-spin asymmetries A\ for electrons and A\ 
for positrons at RHIC, for the various sets of polarized parton distribution functions shown in 
Fig. m Only leptons produced by boson exchange are considered here. 

It is worth pointing out that in Ref. p] spin asymmetries for the same sets of parton 
distributions as in Fig. |6]were shown, but at LO. Our NLO results turn out to be very close to 
the LO ones of [2]. This is easily understood because the bulk of the NLO corrections in the qq' 
channel is the same in the spin-averaged and the polarized cases, so that the corrections cancel 
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to a high degree. This result was also observed in the study [27J of the W cross section without 
the leptonic decay. We do stress, however, that the individual cross sections in the numerator 
and the denominator of Al receive significant NLO corrections of (9(30%). 

For completeness, we show in Fig. [7] the asymmetries computed by including also leptons 
produced by Z/j boson exchange. As expected, the inclusion of 'background' leptons results in 
a reduction of the asymmetry due to the increase in the unpolarized cross section. Consistently 
with the results presented in Fig. [3l the effect is more noticeable at larger rapidities. 
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Figure 7: Same as in Fig. but including also the contribution from Z/7 exchange. 

Figure [8] shows the spin-averaged and single-spin cross sections for electron production, as 
functions of pxi, integrated over \rii\ < 1. As expected, a peak is found around ~ Mw /2 in 
both cases. For the single-polarized cross section, we are again using various sets of polarized 
parton distributions. One can see that the dependence on the polarized distribution functions 
is apparent in the magnitude of the cross section, but hardly in the shape of the transverse 
momentum distribution. The latter is mainly determined by the properties of the gauge boson 
(like its mass and width) and general features of QCD radiation and, therefore, is very similar 
for both unpolarized and polarized cross-sections. In other words, integration over a significant 
region of rapidity washes out many of the features that we found in Fig. O for the polarized 
parton distributions. While this situation may possibly be improved by integrating over non- 
central regions of lepton rapidity, our studies overall confirm our expectation that in order to 
improve our knowledge of the polarized anti-quark distributions it is preferable to study the 
lepton rapidity dependence of the asymmetry instead of its transverse momentum one. 

Apart from providing a good correlation between momentum fractions and lepton rapidity, 
another advantage of rapidity distributions integrated over a wide region of transverse momen- 
tum is that these are insensitive to soft-gluon and resummation effects. As discussed in the 
Introduction, the distribution near pj; = Mvi//2 is sensitive to soft-gluon effects. Proper inclu- 
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Figure 8: Left: transverse momentum dependence of the unpolarized cross section for electron 
production at RHIC. Right: same for the single-spin cross section. Here the lines follow the 
same pattern as in Fig. O 

sion of these effects in tlie theoretical description plays a role, for example, in determinations 
of the W mass from the lepton pj-^ distributions at the Tevatron [21], or in transverse-spin 
studies in W boson production at RHIC [32]. The RHICBOS code [T3|[H] includes an all-order 
resummation of the next-to-leading logarithms in the W transverse momentum qx, which are 
relevant near pti = My//2. However, once one integrates over a sufficiently large region of 
lepton transverse momentum, these logarithms turn into finite higher-order (beyond NLO) cor- 
rections to the transverse-momentum distribution, and neither is their resummation necessary 
nor is it guaranteed to provide an improved theoretical description, as there can and will be 
many other corrections of similar size that are not taken into account. From the point of view of 
extracting polarized parton distribution functions, it therefore seems advisable to us to focus on 
observables integrated over the lepton's transverse momentum, because these are insensitive to 
soft-gluon effects, and to use a plain NLO calculation. Thanks to the fact that the W mass sets 
a very large scale so that the strong coupling is small, and because quark anti-quark annihila- 
tion is the dominant partonic channel, any "partial" beyond-NLO and non-perturbative effects 
remaining from gT-resummation are expected to be relatively small for such observables. Com- 
paring the total W^{W~) spin-averaged cross sections for pt^i > 25 GeV, |?7;| < 1 and setting 
the scales to fJ^F = fi'R = we find 75.4 pb (17.7 pb) with our NLO code, while RHICBOS 
give^ 80.5 pb (18.8 pb). The difference will be in part due to different choices for electroweak 
parameters and parton distributions, but also due to the additional effects in RHICBOS just 
described. Somewhat larger differences between the codes occur if one considers the differential 
cross section at high lepton rapidity and, of course, as a function of transverse momentum at 
Pt,i ~ Mw/2. We have checked that these differences do not, however, significantly affect any 
of the previously performed sensitivity studies by the RHIC experiments. 

■^We thank B. Surrow for providing these numbers. 
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4 Toy Analysis of W Spin Asymmetries in Terms of Po- 
larized Parton Distributions 



While the results presented in the previous section indicate a strong sensitivity of the single- 
spin asymmetries to the polarized light quark and anti-quark distributions, it is quite difficult 
to quantify from them the impact future RHIC measurements will have. In order to investigate 
this, we perform a more detailed analysis. Our strategy is to "simulate" a set of RHIC data 
under hopefully realistic conditions, and to add this set to the data sets included in the published 
DSSV |2l[3] global analysis. A new fit is then performed, for which the simulated data set is 
included, and the impact of future W data from RHIC is gauged from the improvement found 
in this fit for the extracted polarized distributions. As discussed in the Introduction, knowledge 
about the sea quark and anti-quark polarizations in the nucleon so far comes entirely from the 
SIDIS spin asymmetries. In Fig. |9] we recall the DSSV results for the polarized anti-quark 
distributions, including their respective uncertainty bands, which were obtained in DSSV by 
performing a Lagrangian multiplier analysis of the truncated moments ^^^^ dx Au{x, = 
10 GeV^) and J^^^^dx Ad{x^Q'^ = lOGeV^), and allowing modifications of 2% in the total 
of the fit. 



0.04 



0.02 



-0.02 - 




Figure 9: Au{x, Q'^ = 10 GeV^) and Ad{x, = 10 GeV^) as obtained in the DSSV analysis /l 
[^. The bands correspond to changes of 2% in the total of the fit, as discussed in the DSSV 
paper. 

In order to generate a simulated "RHIC data set" we proceed as follows: we compute 
the NLO single-longitudinal spin asymmetries A'l and using the central DSSV set of 
polarized parton distributions □. We then randomly shift the calculated asymmetries, assuming 
a Gaussian distribution of their statistical uncertainties. The latter are estimated using the 
usual formula 6Al = l/{P\fCa). We assume a polarization of P = 60%. For the integrated 
luminosity we consider two values: C = 200 pb~^ and C = 800 pb^^. Concerning rapidity 
coverage, we focus ffist on the present coverage for the Phenix {\rii\ < 0.35) and STAR (|?7;| < 1) 

^Asymmetries are evaluated including both W and Z/7 boson exchange contributions 



14 



experiments. Later, we also investigate the impact of extended rapidity coverage as given by 
\t]i\ < 0.35 and 1 < \r]i\ < 2 for Phenix0 and \r]i\ < 2 for STAR . The generated pseudo-data 
points are shown in Fig. [10] for the two extreme scenarios (smaller luminosity and rapidity 
coverage vs large luminosity and rapidity coverage), along with the predictions according to 
the DSSV set of polarized parton distributions. We have assumed bin sizes in rapidity of 
Arji = 0.33, so that there are six data points in \t]i\ < 1. 




Figure 10: Simulated data generated for e~ (left) and (right) single-spin asymmetries at 
STAR and Phenix for two possible extreme scenarios: integrated luminosity C = 200 pb~^ 
and present RHIC rapidity coverage (upper row), and integrated luminosity C = 800 pb~^ and 
upgraded rapidity coverage as described in the text (lower row). The "simulated data points" 
have been estimated by performing a NLO calculation with the DSSV set of polarized parton 
distributions, followed by a Gaussian shift of the points. The solid lines represent the expectation 
from the DSSV set. 

In order to perform the actual fit, we produce the pre-calculated grids as described in 
Sec. [21 The choice of bins Arji = 0.33 means that we can use the same grids for both RHIC 
experiments, because the rapidity range covered for Phenix corresponds in good approximation 
to two bins with Arji = 0.33. In our first fit, we include in the DSSV analysis the pseudo-data 
generated with the lower luminosity and the rapidity coverage presently available at RHIC. 

^The upgrade planned for Phenix will actually extend the rapidity coverage to 1.2 < |?7;| < 2.4 
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The outcome of this global fit is shown for the polarized anti-quark distributions in Fig. [TTl 
including their resulting Ax^/x^ = 2% uncertainties determined in the same way as in the 
published DSSV analysis. By comparing to Fig. [9l one observes that there is little modification 
of the actual distributions, as compared to the original DSSV ones, but a clear reduction in 
their uncertainty bands. This effect turns out to be very noticeable for x > 0.1, as expected 
considering the rapidity coverage of the pseudo-data added to the global fit. The decrease in 
the uncertainty band is also more noticeable in case of Au, confirming the larger sensitivity 
of e~ asymmetries to this distribution. At values of x ~ 0.01 there is almost no change in 
the distributions and their uncertainties, since the single-spin asymmetries at RHIC are rather 
insensitive to such values of x. The lower row of Fig. [12] shows the result of the corresponding 
fit for the larger luminosity/rapidity coverage scenario. Here the impact of RHIC data can be 
observed down to values of x ~ 0.05, thanks to the extended coverage in rji. 
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Figure 11: Result of a global fit performed by including the simulated data generated with the 
lower luminosity and smaller rapidity coverage scenario (upper row in Fig. [TOj) . 
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Figure 12: Result of a global fit performed by including the simulated data generated with the 
larger luminosity and rapidity coverage scenario (lower row in Fig. [77^) . 
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The pseudo-data used in this analysis were generated to be in full agreement, within sta- 
tistical errors, with the expectation from the DSSV set. Since the latter provides an excellent 
description of the available SIDIS data, we have effectively assumed that constraints on the 
polarized parton distributions emerging from SIDIS and from W production at RHIC are in 
agreement. This should be the case, of course, if both are described by factorized perturbative 
QCD at leading-power. From a theoretical point of view, W production provides the more 
reliable source of information, so if any discrepancy between SIDIS and W production were 
found, it would likely point to higher-twist contributions in SIDIS, or ill-understood issues re- 
lated to fragmentation. In this context, it is interesting to ask what the impact of future RHIC 
data would be if all SIDIS data were removed from the global fit. We have performed such an 
analysis for the scenario with larger luminosity and rapidity coverage. The result is shown in 
Fig. [131 We find that for x > 0.07 the simulated W data put a somewhat better constraint 
on Am and Ad than the SIDIS data presently do |^. Toward smaller x, the distributions are 
of course only very loosely determined because the W spin asymmetries are not sensitive to 
this region. All in all, there are very good prospects for a much better determination of the 
polarized anti-quark distributions from RHIC and SIDIS measurements. 
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Figure 13: Same as Fig. [7^ but excluding all SIDIS data from the fit. 

We end by stressing that there are numerous experimental issues (like efficiencies for lepton 
detection, correct background subtraction, other systematic uncertainties, etc.) that have not 
been included in this simple analysis and that would tend to decrease the impact of the real 
data in the global fit. We regard this study as a "proof-of-principle" that shows that RHIC 
W asymmetry data can be straightforwardly included in a global analysis of polarized parton 
distributions. Future, more detailed, studies will need to be carried out. 

^We remind the reader that the new preliminary COMPASS SIDIS data [8] were not yet included in the 
DSSV analysis fSlj^S] and are hence not included in the present study. 
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5 Conclusions 



We have presented a new next-to-leading order calculation of the cross section and longitudinal 
spin asymmetry for the process pp — )■ i'^X at RHIC, through an intermediate electroweak 
gauge boson. The spin-asymmetry is the main probe of the light quark and anti-quark helicity 
distributions at RHIC. We have developed a multi-purpose Monte-Carlo integration program. 
Our code has the advantage that it allows one to directly include forthcoming RHIC data into 
a global analysis of spin-dependent parton densities, using the Mellin technique of [2|[3|[30]. 
Compared to the RHICBOS code our program does not include any soft-gluon qt- 

resummation effects, as we advocate the use of observables at RHIC that are insensitive to such 
effects. In particular, we have emphasized the advantage of the lepton rapidity distribution over 
the transverse- momentum one. 

Our phenomenological results indicate a good sensitivity of the single-longitudinal spin 
asymmetries at RHIC to the light quark and anti-quark helicity distributions. This finding is 
in line with those of previous studies [2| [T^] - [Ti] . Contributions from Z exchange are found to 
be generally non- negligible. As a benchmark application of our program, we have performed 
a toy global analysis of "simulated" RHIC spin asymmetry data along with the present lepton 
scattering and RHIC high-p^ jet and hadron data. We find that RHIC has a great potential 
for providing better constraints on the light quark and anti-quark helicity distributions, in 
particular at moderately large momentum fractions. Once precise data becomes available from 
RHIC, the consistency and interplay with constraints from SIDIS will be particularly interesting 
to investigate. While more refined sensitivity studies will be needed, we regard our findings as 
a very encouraging signal that precise information on the nucleon's polarized light quark and 
anti-quark distributions will become available before too long. Such information will likely offer 
important insights into the inner structure of the nucleon and the dynamics of QCD. 
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